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In this article we develop a theoretical framework to study the hydrodynamic interactions in the 
presence of a non-fiat and no-slip boundary. We calculate the influence of a small amplitude and 
sinusoidal deformations of a boundary wall in the self mobility and the two body hydrodynamic 
interactions for spherical particles. We show that the surface roughness enhances the self mobility 
of a sphere in a way that, for motion in front of a local hump of the surface, the mobility strength 
decreases while it increases for the motion above a local deep of the rough surface. The influence of 
the surface roughness in the two body hydrodynamic interactions is also analyzed numerically. 
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I. INTRODUCTION 

Hydrodynamic interaction of colloidal particles in con- 
fined geometries is an important problem in low Reynolds 
fluid dynamics [T| . Fluid motions confined by one or two 
parallel flat planes are interesting examples with either 
analytical or numerical known solutions 0, H| . The so- 
lution to these problems are involved in soft matter re- 
lated phenomena and also microfluidic experiments ■ 
In soft matter systems, swimming motion in a geomet- 
rically confined environment is a subject of growing in- 
terest 0413 , apparently better understanding of these 
systems requires a good knowledge of the hydrodynamic 
interactions in confined geometries. In microfluidic ap- 
plications, a better control of the processes requires the 
prediction of the hydrodynamic effects due to the walls. 

Confining wall can generate nontrivial effects. For ex- 
ample, in a very simple system composed of a single 
sphere near a wall, the mobility parallel to the wall is 
always larger than the perpendicular mobility [l| which 
has been verified experimentally (Til [l2j ■ As other ex- 
amples for the effects due to the boundaries, we address 
the experiments, showing that microorganisms, e.g. E. 
Coli (Ta |. bull spermatozoa [TjJ, swimming in confined 
geometries are attracted by surfaces. In this article we 
concentrate on the effects due to the roughness of the 
confining walls. We will consider a rough wall that con- 
fines the fluid flow at low Reynolds number, and ask the 
following question, does the long wavelength roughness of 
the wall have any important influence on the one or two 
body hydrodynamic interactions? We use a perturbation 
method and investigate the case of a small amplitude and 
regular roughness on an infinite wall. We show that the 
roughness has important contributions in the hydrody- 
namic interactions. 

Another challenging issue in the low Reynolds, quics- 
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cent fluid dynamics, is the validity of no-slip boundary 
condition. This is certainly important where the nano 
structure of the surface is involved. There are experi- 
mental and theoretical works, investigating the influence 
of the nano-roughness of the boundaries a nd j ustify the 
validity of the so-slip boundary conditions jl5l. [l6j. Here 
we would like to stress the fact that, justification of the 
no-slip boundary condition, necessarily needs analytical 
results for rough surfaces taking into account the no-slip 
boundary condition. 

The rest of this article is organized as follows: In Sec. 
II, we present a short review on Stokes flow and introduce 
the hydrodynamic interactions. Section III is devoted to 
the hydrodynamic effects of a rigid and flat wall. The 
effects of a general roughness is presented in sec. IV. An 
example of a sinusoidal roughness is considered in Sec. 
V. Concluding remarks are presented in sec. VI. 

II. STOKES FLOW AND HYDRODYNAMIC 
INTERACTIONS 

To study the fluid motion for colloidal particles sus- 
pended in a fluid medium we define the Reynolds num- 
ber as the ratio between the characteristic transport time 
scale due to diffusion and the convection time over a 
length L. Denoting the fluid density by p, viscosity by 77, 
typical velocity by U and also the linear size of the parti- 
cles by a, we see that the Reynolds number is Re = I^2£. 
In a wide variety of phenomena occurring in the motion 
of micron scale particles, the Reynolds number is very 
low. For these phenomena we can consider the limit of 
zero Reynolds number. At zero Reynolds number the 
fluid dynamics is expressed by Stokes equation. Denoting 
the fluid velocity and pressure fields by u(x) and P(x), 
the Stokes and continuity equations for an incompressible 
flow can be written as: 

77V 2 u(x)-VP(x) = f e2;t , V-u(x) = 0, (1) 

where f ext ; denotes the density of external body force 
acting on the fluid. The fluid velocity field is subject to 
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FIG. 1: (color online). Two spherical particle is moving in a 
low Reynolds flow. The fluid flow is bounded by a rough, rigid 
and no-slip boundary. Velocities and hydrodynamic forces 
acting on particles are shown in this picture. 

no-slip boundary condition on the surface of solid bound- 
aries. We would like to consider the motion of N colloidal 
particles suspended is a low Reynolds flow (see Fig. [T]). 
Recalling the fact that the governing equations are lin- 
ear with respect to the velocity profile, and applying the 
no-slip boundary condition on the surface of rigid par- 
ticles, we can eliminate the fluid velocity and obtain a 
set of linear equations relating the particle velocities to 
the hydrodynamic forces acting on them. Denoting the 
velocity of m'th particle by v m and the corresponding 
hydrodynamic force by f m , we can write the following 
equation: 

N 3 

< m>v" • 'T- ( 2 ) 

n=l j=l 

Coefficients -D™" are know as the elements of hydrody- 
namic interaction. Quite similar to the Onsager relations 
in thermodynamics and based on general symmetry ar- 
guments and microscopic reversibility, the hydrodynamic 
tensor £>J" n is symmetric with respect to its up or down 
indices [I3,[ll] ■ For special case of 2 particles (denoted by 
a and /?), and for convenience we define the self mobility 
tensor as the response of a particle to the force acting on 
it by fiij = Df 3 a . And also the hydrodynamic response of 
a'th particle to the force acting on /3'th particle is shown 
by: Mjj = D"? = D^ a . Different components of the hy- 
drodynamic interaction tensor depend on the size and the 
relative configurations of the particles. Calculating the 
hydrodynamic interaction tensor, is an important prob- 
lems in low Reynolds hydrodynamics. Here we present an 
analytic method that is useful for very small and spher- 
ical particles. A very small sphere can be considered as 
a point force, a source term in the Stokes equation. De- 
noting the point force by f(x) = b <5(x — xA the velocity 
field of the point force, can be written as [19L [20l| : 

3 

u *( x ) = ^2 G ij( x > x o) x bj, (3) 

i=i 

and correspondingly the associated pressure field is given 
by: P(x) = II(x, xq) • b. Here Gy(x, Xq) is the Green's 



function of the Stokes differential equation. 

Solving the Stokes equations and having the associated 
Green's function for any required geometry, we can use 
Faxen's theorem for spherical objects with radius a to ex- 
press the mobility tensor in terms of the Green's function 
of Stokes equation P, [2l| : 

Dff (x^) = (i + y v»„)(i + ^vi^G^yy 

(4) 

For very small spheres with radius a smaller than the 
typical distance between spheres, then Dfj 1 (x Q , x* 9 ) = 
G^(x Q ,x' 3 ) up to second order in a [22l |. 

In the rest of this article, we use the above method and 
presents the results for unconfined ([/), confined with a 
fiat wall (F) and confined with a rough wall (R). The 
effects due to roughness is considered for a very small 
amplitude roughness. 

III. HYDRODYNAMICS NEAR A FLAT AND 
NO-SLIP WALL 

We first consider the case of a fluid flow that is bounded 
by a rigid, flat (F) and no-slip plane located at z = 0. 
The solution to the problem of a point force in the pres- 
ence of a flat wall is given by Blake who has used 
the image method to construct the solutions. The flat 
wall Green's function for the upper half space (z > 0), is 
given by: 

G£(x,x ) = Gg(x,xo)-Gg(x,xjr) 

+ C^(x, x™) - Gy D (x, x™), (5) 

where the Green's function for an unbounded (U) fluid 
flow GVj, that is called stokeslet is given by: 

here r = x — x , and x m = x — 2 (z ■ x ) z is the image 
position of the point force with respect to the flat wall. 
Defining a new vector R = x — x m , we can express the 
potential dipole field G^, and stokeslet dipole G^ D as: 

aS^)-^*i-«,.)^(§), (7) 



G| D (x, x ;"') = 22b(l - 2«„) M -Gg(x, x <™), (8) 

where z — z.x. The pressure field for the flow that is 
bounded by a flat wall is given by: 

n F (x,x ) = n l, (x,x ) + n c/ (x,x m ) 

- 2z 07 J r (^n c/ (x,xn), (9) 



3 



where the pressure field associated to a point force in an 
unbounded space is given by: 



n^(x,x ) = — 



2 (x - x ) 



x 



(10) 



Having in hand all the above information, we can calcu- 
late the hydrodynamic interactions in the presence of a 
flat no-slip wall. 

As an example and using the above formalism we 
present the results for the self mobility of a very small 
sphere moving above a flat and no-slip wall. Denoting 
the perpendicular distance between the sphere and wall 
by H, and assuming that H is greater than the sphere 
radius a, the components of self mobility tensor are |i1.[2lj| 



F 



Mo(l 
Mo(l 



16 H ( H> h 
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where fio = l/(67rr?a) is the self mobility of a spherical 
particle moving in an unbounded fluid. As an important 
and non-trivial result, the above equations show that, the 
self mobility of the sphere in the direction parallel to the 
wall is always larger than the perpendicular direction [l| . 
This effect has been verified experimentally 11 1. Calcu- 
lations show that, up to O(-g-) 2 , other components of the 
self mobility tensor are zero. 



IV. EFFECTS DUE TO A ROUGH AND 
NO-SLIP WALL 

The aim of this article is to express analytical expres- 
sions for hydrodynamic interactions in a semi infinite 
flow, bounded by a rough (R), rigid and no-slip wall. 
Fig. Q] shows the schematic view of a rough plane that 
bounds the fluid flow. Two spherical particles located at 
positions x Q and x' 3 are moving in the fluid. The position 
vector for the points on the bounding wall is considered 
as: x s = (x, y, h(x, y)). This kind of parametrization al- 
lows us to obtain the results of a flat wall by considering 
the limiting case of h(x, y) — 0. For later use we define 
x" = (x, y, 0), that is denoting the location of the points 
of a flat wall located at z — 0. 

Applying the no-slip boundary condition on the rough 
surface (R), we expect that the influence of the rough- 
ness will produce non-trivial effects. To study the hy- 
drodynamic effects of a rigid and rough wall with small 
amplitude roughness, we can construct a perturbation 
expansion. Introducing a small dimcnsionless parameter 
e = ho/z, where ho is the typical amplitude for height 
fluctuations and z measures the distance of particles from 
the wall, we expand all quantities in powers of e. As ex- 
plained before, hydrodynamic interaction between small 
spherical particles can be obtained by solving the velocity 
field of a point force for the required geometry. Here the 
solutions to the Stokes-Green's equation should be ob- 
tained by applying the required boundary condition. The 



velocity field of a point force satisfies the Stokes equation 
(Eq. [I} and is subject to to the following boundary con- 
dition: 



u(x s ) = 0. 



112) 



Expanding the corresponding velocity and pressure field 
of a point force in the presence of a rough wall in powers 
of e, we can write: 

u fl (x) = u(°)(x) + u«(x) + u( 2 )(x) + 0(£ 3 ), 
P fl (x) = P(°)(x) + P«(x) + P( 2 )(x) + 0(£ 3 ),(13) 

The zeroth order term are the velocity and pressure field 
of point force in the presence of a flat(P) and no-slip wall. 
In this case we will have: 



u (0) (x) = u F (x) , (x) = P F (x) 



(14) 



Higher order corrections due to the wall roughness, can 
be obtained by noting that different order of the velocity 
and pressure fields are satisfied the following differential 
equation: 

77V V' 1 ' (x) - VP (n) (x) = 0, V ■ uW (x) = 0, (15) 

where the boundary conditions for the first and second 
orders are given explicitly by: 



u (2) (*°) 



^dJ)(y°) ft 2 (x°)d 2 u<°)(x°) 



dz 



d 2 : 



(16) 



As one can see the velocity field at order (n), is related 
to the local value of the derivatives of the velocity filed 
at order of (n — 1), in the position of flat wall. Using 
the well known integral representation for the flow field 
of Stokes equation, we can write the following integral 
solution [l9j ): 



U ? } W = ^ / uj"V)T 4 , fc (x,x')dS fc (x') 



(17) 



where the integration is carried out on a closed surface 
composed of an infinite plane at z = 0. The surface 
is closed at the upper half space and dS(x') is the area 
element on this surface in the inward direction (here is 
I). The stress tensor Tjjj(x,x') is given by: 



T ijk (x, x') = -6nS ik Uf (x, x') + Stttj^-G^ (x, x') 



8vrr?— G^-(x,x'). 



(18) 



Now we can expand the hydrodynamic interaction tensor 
moving in a fluid bounded by a rough wall, in terms of 
small parameter s as: 



T> R = D<°) + D« + D< 2 > + 0(e 3 ), 



(19) 
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where the zeroth order is the results of a sphere moving 
in a medium bounded by a flat (F) and no-slip wall: 



1^ = Df 3 



G 



i j ' 



(20) 



Using Eqs. (J3|), (@|, (fTT|) and considering the explicit 



form of the stress tensor TV^x^Xq), we can arrive at 
the following expressions for the first and second order 
corrections to the mobility tensor due to the roughness 
of the wall: 



StttjxD^x^) = -i- |^ )(x Q ,xO) J D'(°)(x /? ,xO)/ l (xO)d 2 xO 



(21) 



256tt 2 77 x L>( 2) (x Q ,x fi 



1 

2^ 



^(°)(x Q ,x2)^(°)(x,,x'°)I?"( )(x°,x^)Mx°)Mx^)d 2 x°d 2 x^ 
- / (^(x Q ,x2)Z)(°)(x /j ,x2) + (a^P))h 2 W° s , 



(22) 
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FIG. 2: (color online). A spherical particle is moving in a con- 
fined fluid flow. The roughness of the boundary is assumed to 
be a sinusoidal pattern with wavelength A in the x direction. 
The boundary is assumed to be invariant under translation 
along the y direction. Two parameters, H the average verti- 
cal distance from the wall, and w the vertical deviation from a 
local hump of the surface roughness, characterize the sphere's 
position. 



where 



D'(°)(x^x°) = ^Wx) 



^" (0) (x°,x'°) = 



d 2 



dzdz 
d 2 



7 L> (0) (x,x') l(x,x')-Kx°,x'°) 



£»^(x^x°) = _£>(°)( X/3 ,x) U x o, 
oz^ 



and = Tl\p6 i=z + D^' ' 5^ z . As one can see, the 
mobility tensor is symmetric: D(x Q ,x ( g) = D(x.p,x. a ). 

In the following section we choose a very special form 
of the wall deformation pattern and investigate the hy- 
drodynamic interactions. 
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V. SINUSOIDAL ROUGHNESS 

As a special example we assume that the wall rough- 
ness is a sinusoidal deformations of a flat wall located 
at z = 0. The height profile of the rough plane, in the 
Monge representation, is given by the following equation: 



z = h(x, y) = ho cos(qx 



(24) 



where ho is the amplitude of the roughness, and the wave 
vector of the roughness along the x axis is denoted by q. 
This special type of roughness allows us to study any gen- 
eral roughness that is invariant under translation along 
y direction. 

As a special feature of the hydrodynamic interaction 
we concentrate on the self mobility of a spherical parti- 
cle moving in a medium bounded by a rough and no-slip 
wall. Fig. [3J shows the schematic view of a sphere with 
size a, moving near a rough wall. Perpendicular pro- 
jection of the sphere position into the z — plane, is 
deviated from the wall hump by a distance w. Denoting 
the wave length of the wall roughness by A = ^ , and the 
average height of the particle on the wall by H , we define 
a dimensionlcss parameter by 7 = 27Ty. Now following 
the method developed in the preceding sections, the first 
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and the second order corrections are given by: 
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8ttt]H H 2 I 32 

-^o(27) 
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fcos(20)[^ 7 3 if 3 (27) 
|7 5 A"i(27) 



- 7 ^ 1 (2 7 )-9 7 t 



+36cos 2 (0)e~ 7 (l -7) 



5 J (2s) 



ds 



(s 2 + 7 2 ) 5 
/ 2 ( 7 /2)^ 2 (7/2) 
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ryy 
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8iri]H H 2 132 32 



cos(20) 7 3 ^ 3 (27) 



-36 cos 2 (</>)e~ 



1 



8irr]H H 2 I 8 
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7 / 2 (7/2)A- 2 ( 7 /2)}, 

+ | cos (2^)[_ 7 4 jft r 4 (2 7 ) 
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7 5 A' 5 (2 7 )]}, 



(26) 



Where = measures the deviation of the perpen- 
dicular projection of the sphere from a local hump on 
the surface. Here J n (l) , Jn(l) and K n {"f) are Bessel's 
function of the order n [231 ] . 

We can investigate the effects of the roughness in two 
extreme limits of long or short wavelength deformations. 
For the case of very short wave length deformations, A <C 
H (7 3> 1) one can see that the roughness has no net 
effects on the self mobility of a spherical particle moving 
very far from the wall (up to first order of e). In this case 
the mobility tensor is effectively given by the mobility 
tensor of a sphere located very far from a flat wall. The 
explicit form of the mobility tensor for this case is given 
by Eq. (fTTj) . The back flows, scattered from the humps 
and deeps of the wall have cancelled out the effects of 
each other and the overall averaged back flow looks like 
a back flow from a flat wall. 

In the case of a very long wavelength roughness where 
A > (7 C 1), we proceed and obtain the second order 
corrections (in terms of e) to the self mobility compo- 
nents. Now we can expand this result around small 7, to 
reach the following expressions for the different elements 
of the self mobility tensor: 



r-xx 



u R 

ryy 



f^zz 
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^X2 
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5cos(20)) +3cos 2 (0) 



/4 + Mo§(^)(l + cos(20)). 



(27) 



As one can see, the effects due to the roughness, en- 
hances the self mobility tensor of a sphere in asymmetric 
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FIG. 3: (color online). Different components of self mobility 
tensor 5 fin = (/if| — ^fi)/^fi for a spherical particle moving 
adjacent to a non-smooth and no-slip wall are plotted as a 
functions of the particle's separation from the wall. Here the 
wave vector for surface roughness lies along the x direction. 
All three components of the self mobility tensor are plotted 
for two special cases of motion above a hump or deep of the 
surface deformations. As one can see for = (sphere move 
over a local hump) , the roughness decreases the parallel com- 
ponents of self mobility while for cj> — t (motion over a local 
deep), they increase. Parameter values for these graphs are: 
ho/H = 0.1, a/H = 0.1. 



way with respect to the in plane (x and y) directions. 
We define the asymmetric parameter as the difference 
between the mobilities in the x and y directions: 



= V>xx - Vyy = A'O 



21 , ah 0\ f±\ 



(28) 



The sign of the asymmetric parameter depends on the 
local position of the sphere. In terms of the back flow 
scattered from the wall and in the limit of long wave 
length deformations, it is expected that the nearest hump 
or deep will have the dominant contributions on the parti- 
cle motion. Over a local hump (0 < 4> < 7r/2), \i r x > [i R y , 
while over a local deep (tt/2 < (f> < 3ir/2), fj, R x < [i R y . 
Interestingly for a special points of <fi = 7r/2,37r/2, the 
mobility tensor is symmetric. 

To analyze the self mobility for intermediate 7 {H ~ 
A), we have presented numerical plots in Fig. [3J As one 
can see for H ~ A, the change in the mobility of a sphere 
mediated by roughness changes its sign for motion above 
a peak or above a valley. 

To demonstrate the effects of roughness on the two- 
body interactions, we have presented in Fig. 0] and Fig. 
[SJ numerical results for diagonal and off diagonal ele- 
ments of hydrodynamic interactions. In these examples 
we have studied the hydrodynamic interactions for some 
different cases. In Fig. 2] we have plotted M xx (r) (up) 
and M yy {r) (down). In Fig. @] right part (up and down), 
we have plotted the interactions of two spheres that have 
same vertical distance H from the wall. As one can see, 
the roughness always decreases the strength of M xx (r) 
and M yy (r). There is a very weak periodicity with the 
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FIG. 4: (color online). Different components of the hydrodynamic interaction tensor 8Mu = 8irr)H(M^ — M^), for two 
spherical particles moving parallel to the wall, are plotted as a functions of the particle's separation. Here the wave vector for 
surface roughness lies along the x direction. At right graphs (up and down), the interaction is plotted for the case where the 
position of the first sphere is fixed at (f> = and we changed the position of the second sphere. There is a periodicity with the 
wavelength of the roughness. As one can see in the case of left graphs, for <f> = (sphere move over a local hump), the roughness 
decreases the interaction strength while for <f) = n (motion over a local deep), it increases. Parameter values for these graphs 
are: \/H = 10, ho/H = 0.1, a/H = 0.1. 
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FIG. 5: (color online). Off diagonal components of the hydro- 
dynamic interaction tensor SMy — 8Tir)H (M^ — My) , for two 
spherical particles, are plotted as a functions of the particle's 
separation. Here the wave vector for surface roughness lies 
along the x direction. Parameter values for these graphs are: 
\/H = 10, h /H = 0.1, a/H = 0.1. 



wavelength of the roughness. In the left part of Fig. @] 
(up and down) , we have plotted the hydrodynamic inter- 



action as a function of vertical separation of two spherical 
particles. The results for this case, depend on the local 
position of the spheres with respect to the surface rough- 
ness. As one can see for <fi ~ (spheres move over a local 
hump), the roughness decrease the interaction strength 
while for <f> = tt (motion over a local deep), it increases. 
As another example, some of off diagonal components of 
the two body hydrodynamic interactions are plotted in 
Fig. [5] There is a very weak periodicity with wavelength 
A, that is not clearly seen in the scale of this graph. 



VI. CONCLUDING REMARKS 

In this article we have considered the influence of a 
rough, rigid and no-slip boundary wall on the hydrody- 
namic interactions of spherical particles. We have stud- 
ied a regular sinusoidal roughness pattern with very small 
amplitude roughness on a flat plane. For simplicity we 
have studied a simple wave with a single wave vector 
along x direction. Taking into account the wall effects 
by applying the no-slip boundary condition by standard 
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perturbation technique, we have calculated one and two- 
body hydrodynamic effects. For a single and small ra- 
dius sphere moving in the presence of a rough wall, we 
show that, the different elements of self mobility ten- 
sor changes in asymmetric way. Motion along the wave 
vector of surface roughness is different from the other in- 
plane direction. When a spherical colloid suspended near 
a rough and no-slip wall, the hydrodynamic drag force 
depends on the local position of the sphere. Roughness 
will produce different contributions for motion on a local 
hump or a local deep of the wall. This kind of behav- 
ior in two-body hydrodynamic interaction is also seen by 
numerical investigations of different components of the 
hydrodynamic interactions. 

We note that in the current formulation of the problem, 
we have made some approximations which need to be 
dealt with carefully. Caution is needed in applying the 
results since there are many length scales in the problem: 
a, ho, A and H . First, the Faxen's formula has allowed us 
to treat the sphere's size in a series expansion in powers 
of ei = (a/H). Second approximation is related to the 
slowly varying roughness of the wall and consequently 
to a series expansion of the results in powers of ti = 
(ho/H). The series expansion for a typical component 
of the hydrodynamic interactions, for example the self 
mobility, have the following structure: 

A M /Mo=eiX { [fu(r/)e 2 + fi2(r/)4 + 0{4)] 

+ [/ 8 i(7)ea + f22h)e 2 2 + 0{efj] e l 
+ 0{e\)), (29) 



where we have already defined 7 = 2tt(H/X). The results 
are valid for any 7, however one should note that, the 
convergence of the above series expansion in the limit 
of very small 7 (7 1), is constrained to the condition 
(ho /A) <C 1. By small roughness assumption, one expects 
that this criterion is satisfied as well. 

In conclusion, we have developed a systematic way 
to evaluate the perturbation effects of a rough wall in 
the hydrodynamical properties of small spherical parti- 
cles. The results of current work can be used in many 
directions related to the colloidal problems in confined 
flows, where the roughness is an ignorable characteristics 
of most walls. The influence of wall roughness on the 
thermal diffusion of colloidal particles is an interesting 
issue that we are considering. Inspired by the ensem- 
ble of low Reynolds swimmers, we are also investigating 
the motion of low Reynolds self propellers adjacent to a 
rough wall. 
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